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This work is concerned with phrasing the concepts of fault-tolerant quantum computation within the frame- 
work of disordered systems, Bernoulli site percolation in particular. We show how the so-called "threshold 
theorems" on the possibility of fault-tolerant quantum computation with constant error rate can be cast as a 
renormalization (coarse-graining) of the site percolation process describing the occurrence of errors during 
computation. We also use percolation techniques to derive a trade-off between the complexity overhead of the 
fault-tolerant circuit and the threshold error rate. 
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Many researchers' confidence in the eventual experimental 
realization of reliable large-scale quantum computation rests 
upon a number of results known collectively as the "thresh- 
old theorem" (see, e.g., d 011 H). The recurring motif of 
these theorems is that, under certain reasonable assumptions 
on errors in the computer, fault-tolerant quantum computation 
is possible provided that the error rate rj does not exceed some 
threshold value rjc- 

Among the papers dealing with threshold theorems, the 
work of Aharonov and Ben-Or 01 is especially noteworthy: it 
contains rigorous constructive proofs of the possibility of arbi- 
trarily reliable sub-threshold quantum computation for a wide 
variety of noise models, including local probabilistic noise 
(i.e., when each gate in the computer suffers an error with 
probability 77 and functions correctly with probability 1 — 77, in- 
dependently of all other gates both in space and in time), noise 
with exponentially decaying space-time correlations, and gen- 
eral noise (i.e., one not described by an a priori probabilistic 
model). The threshold theorem is also specialized to quantum 
computation on the d-dimensional hypercubic lattice Z'* with 
a restriction that only nearest-neighbor qubits can interact di- 
rectly. 

The strategy of Aharonov and Ben-Or has a lot in common 
with the work of Gacs jsl on fault-tolerant classical computa- 
tion in cellular automata. In particular, at the core of their 
proof lies the idea of iterated simulation of one unreliable 
computer by another, with (quantum) error correction imple- 
mented on all levels of iteration. It is then shown that, pro- 
vided that the eiTor rate rj (say, per quantum gate) is below a 
certain threshold rjc, each iteration will reduce the effective 
error rate (rate at which errors occur in the encoded infor- 
mation). One of the goals of this paper is to provide an al- 
ternative interpretation of the Aharonov-Ben-Or method for 
local probabilistic noise in terms of a simple disorder model, 
namely Bernoulli site percolation In particular, we will 
exhibit a close relation between the recursive simulation tech- 
nique and a renormalization (coarse-graining) of the site per- 
colation process, and then use this relation to derive a trade- 
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off between the threshold error rate rjc and the complexity- 
theoretic overhead required to implement the computation 
fault-tolerantly, i.e., the minimum number of iterations needed 
to bring the effective error rate down to some desired level. 

Let us sketch very briefly the key ideas behind iterated sim- 
ulation ||4i]. The basic ingredient is a quantum error-correcting 
code Q], i.e., an isometric embedding of a Hilbert space Hk of 
k qubits as a fc-qubit subspace IC of some Hilbert space Hm of 
TO qubits, with m > k [this is refeiTed to as a quantum (to, k)- 
code]. Aharonov and Ben-Or 0] use quantum {m, l)-codes. 
Errors are modeled by quantum operations |7], i.e., mappings 
of the form T(p) = EjpE]^, where p is a state (density 
matrix) on TLm, and the operators Ej : Hm Ti-m (called the 
&flM.soperatorsofT) satisfy the constraint E'jiJj = l. A 
quantum {ni, l)-code is said to correct d errors (with d < m) 
if there exists a quantum operation R (the recovery opera- 
tion), such that for any density matrix supported on the code 
space IC and for any operation T whose Kraus operators are 
tensor products of at most d nontrivial single-qubit operators 
acting on the qubit components of K, and identity operators 
on the rest of the qubits, we have R o T(p) = p. If this is 
the case, we say that the code Hk ^ IC Hm is a quantum 
{■m,k, d)-code. 

When quantum error-coiTecting codes are used to protect 
information inside a quantum computer, each qubit is encoded 
separately. Thus, if we use a quantum {m, l)-code and there 
are n qubits, then the encoded state of the computer's register 
is a state of 0{mn) qubits (the constant hidden in the "big-oh" 
notation reflects the ancillae added to each ?7i-qubit block in 
order to implement the recovery operations). We assume that 
all encoding, decoding, recovery, and computation operations 
are implemented using quantum gates from a suitable univer- 
sal set 1 8]. Aharonov and Ben-Or use concatenated codes. 
That is, each qubit is encoded in a block of to qubits, each of 
which is in turn encoded in a block of to qubits, and so on. If 
we do r levels of concatenation, we end up with n blocks of 
rrf qubits each. At the end of computation decoding proceeds 
hierarchically, starting with the highest (coarse-grained) level 
and ending with the lowest (fine-grained) level. 

The computation itself is also defined hierarchically. Con- 
sider first the case of r = 1 (one level of concatenation). Each 
gate J7 G 5 in the original circuit now corresponds to a cer- 
tain composition of gates from Q, referred to as the procedure 
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of U . Each time step of the original (unencoded) computa- 
tion now corresponds to a working period, consisting of two 
stages: the recovery operation appUed on each m-qubit block, 
followed by parallel application of the procedures of all the 
gates used in the original computation in this particular time 
step. Denoting the original quantum circuit by Co and the en- 
coded circuit by Ci, we have a mapping M : Cq ^ Ci. This 
is, essentially, a simulation of the original circuit Co by the 
encoded circuit Ci . 

Now imagine the same construction done with the circuit 
Ci, treating each m-qubit block as a unit (a 1-block, in the ter- 
minology of Aharonov and Ben-Or |4]), and each procedure 
as a single gate. This gives a circuit C2 = Af (A/(Co)) = 
Af^(Co). After r levels of concatenation we will end up with 
a circuit Cr = AI^{Cq) that simulates the circuit C^^i that 
simulates the circuit Cr-2, and so on. The coarse-grained 
computation in Cr takes place on n r-blocks, each of which 
consists of m'' qubits on the fine-grained level. Again, at the 
end of computation decoding takes place recursively, starting 
at the level of r-blocks and so on, until the level of individual 
qubits (0-blocks) is reached. 

Note that the concatenated circuits thus constructed are es- 
sentially self-similar. That is, if we start with Cr-s and do s 
levels of concatenation, then the resulting circuit Cr will "look 
like" the circuit Cs = Af'*(Co) if each block of rn" qubits in 
Cr is treated as a unit (an s-block 0]). 

Finally, we have to safeguard the computation against the 
propagation of errors. The idea is to design all encod- 
ing/decoding operations and all procedures in such a way that 
they are fault-tolerant. Informally, this means that a small 
number of localized errors during a procedure will not affect 
too many qubits at the end of the procedure, so that the recov- 
ery operation which is applied during the next working period 
will be successful with high probability |!4]. More precisely, 
given a quantum (to, 1, (i)-code, we say that it has spread s if 
a single error anywhere during a procedure results in at most 
s faulty qubits in each block at the end of the procedure. Note 
that it follows from self-similarity that this definition of the 
spread of the code makese sense at any level in the concate- 
nation hierarchy. Now if the code under scrutiny is a quantum 
(to, 1, d)-code, then we require that s < d, so that at least 
one error can be tolerated in each procedure. Aharonov and 
Ben-Or |4] call such codes quantum computation codes. 

In order to visualize the occurrence of errors during compu- 
tation, we will associate an interaction graph to the quantum 
circuit under consideration. Given a circuit C, we will con- 
struct its interaction graph Tc in two steps. First, we will re- 
place each quantum gate in C, including the identity gates, by 
a single vertex. This intermediate combinatorial object will 
in general be a multigraph, i.e., there can be multiple edges 
joining a given pair of vertices. The second step is to col- 
lapse each such bundle of multiple edges into a single edge. 
This construction is illustrated for a simple quantum circuit in 
Fig-E In precise terms, the interaction graph corresponding 
to the quantum circuit C is given by Vc = {Vc, Ec), where 
the vertex set Vc is the set of all quantum gates, including the 
identity gates, that are used in the computation, and the edge 
set Ec consists of pairs {v, v'), v, v' £ Vc, such that there 




FIG. 1 : A simple quantum circuit and the corresponding interaction 
graph. The solid vertices correspond to nontrivial quantum gates, the 
empty vertices correspond to identity gates. Thick edges are those 
that resulted from the collapsing operation on the intermediate multi- 
graph. 

is at least one quantum wire in C connecting the gates corre- 
sponding to V and v' . We will use the shorthand v ^ v' io 
denote the fact that {v,v') € Ec. 

The vertices of Tc can now be thought of as the potential 
locations of errors, so that an error that occurs at some v G 
Vc can propagate to some or all of v' £ Vc for which v <^ 
v' . Assuming that there are no correlations between errors in 
distinct gates, either in space or in time (which is the focus 
of this paper), we lose no generality if we suppose that Tc is 
connected, i.e., there exists a path between any two v, v' S 
Tc. (Otherwise we could represent Tc as a union of maximal 
connected subgraphs, Tc — T^ U . . . U T^, so that if an error 
occurs in some Tq, 1 < i < L, it will not propagate to any 
T-jj with j ^ i.) Another important feature of the graph Tc 
is that, even if it corresponds to an "infinitely large" quantum 
circuit, its vertices have bounded coordination number, i.e., 

max ^ {{v' G Vc I w ^ v'}) < 00, 

v&Vc 

where #(5) denotes cardinality of the set S. Basically, this 
is a consequence of the fact that all quantum gates in any uni- 
versal set act on finitely many qubits. 

As for the error model, we will content ourselves with the 
simplest scenario. Namely, some rj G (0, 1) is picked, and we 
assume that each gate (including identity gates) undergoes an 
error with probability 77 and functions correctly with probabil- 
ity 1—77. This is known as the assumption of local stochastic 
faults 1 4]. Quantum operations that describe errors within this 
model have the form T(p) = (1 — ri)p + riT{p), where T' 
is a quantum operation whose Kraus operators act nontrivially 
only on the qubits participating in a particular quantum gate. 

Under this simple model, the occurrence of errors in 
space and time is described naturally by a process known as 
Bernoulli site percolation |6] on the interaction graph Tc. 
Namely, each vertex w G Vc is occupied with probability rj 
and vacant with probability 1 — rj, independently of all other 
vertices, so that the occupied vertices correspond to the loca- 
tions of errors. Configurations of occupied and vacant vertices 
are given by elements of the sample space ^c ■— {0, l}^*^. 
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where (1) indicates a vacant (occupied) vertex; the proba- 
bihty of an event F C 51^ will be denoted by (F) [note 
the explicit parametrization of all probabilities by the vertex 
occupation density rj]. An edge v ^ v' is called open if both 

V and v' are occupied, and closed otherwise. A set K <Z Vc 
is called a cluster (of connected occupied vertices) if all edges 

V ^ v' with v,v' G K are open, and all edges v ^ v' with 

V ^ K,v' ^ Vc\K are closed. Among other things, percola- 
tion theory is concerned with statistical properties of clusters 
as the occupation density rj is varied. 

By construction, the interaction graph Tc and the corre- 
sponding site percolation process have the following proper- 
ties: (i) Tc is statistically homogeneous, i.e., the probability 
that a vertex is occupied does not depend on its location in 
Fc and (ii) the number of vertices that can be reached from 
a given vertex by paths of length n or less is 0(2^") for some 
c < 1. Then it follows from general arguments of percolation 
theory {S\ that: (1) there exists a number 77* > (the percola- 
tion threshold |9]) such that, in the limit of an infinite number 
of vertices, the probability of a given vertex belonging to an 
infinite cluster of occupied vertices is zero for 77 < ry,, and 
strictly greater than zero for 77 > 77,; and (2) for 77 < 77,, the 
expected cluster size is finite. 

We now show that an estimate of the threshold error rate 
may be derived by means of a renormalization argument com- 
mon in percolation theory (see, e.g.. Chap. 4 of |6] for de- 
tails). The basic idea is the following: an (infinite) graph F 
is "coarse-grained" by means of a rule that replaces suitably 
defined groups of vertices of F with single vertices, prescribes 
how these new vertices are to be connected by edges, and de- 
fines a new occupation probability 77' on the resulting graph 
F'. Typically, if we start with a subcritical (77 < 77*) percola- 
tion proces on F, the net result of iterating the renormalization 
process will be to push the system away from 77* towards the 
"trivial" behavior (77 0). 

Given a quantum circuit Cq, consider the single-layer en- 
coding Ci = Af (Co) and the accompanying quantum com- 
putation code. Let A be the maximum number of locations 
involved in a single procedure in Ci , d the maximum number 
of errors corrected by the code, s the spread of the code, and 
k = \d/ s\. Fix a particular procedure in Ci. Then, if more 
than d qubits are in error at the end of this procedure, the re- 
covery stage of the subsequent procedure will fail. This will 
happen precisely when at least fc + 1 errors occur during the 
procedure. Denote this event by E. Using subadditivity of 
the probability measure P,,(-) and properties of the binomial 
coefficients, we readily obtain the bound 
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Then, for any 77 satisfying the threshold condition rj < i]c = 
we have F,^{E) < rj, i.e., the probability of k or more 
errors during a procedure is smaller than the probability of a 
single error. 

The influence of concatenation on the effective error rate 
can now be understood as follows. Consider the interac- 
tion graph associated with the circuit Ci = M{Co). 



We "renormalize" it by replacing the vertices corresponding 
to each procedure in Ci with a single vertex, and then by 
drawing edges appropriately. From considerations of self- 
similarity we can expect that the resulting renormalized graph 
i?(Fci) is isomorphic to Tco- We then say that a vertex v 
of -R(F(7i ) is occupied if at least fc + 1 errors occur in the 
corresponding procedure in Ci, and vacant otherwise. We 
will denote the occupation density in i?(Fcj) by i?(7/); it 
is related to the "original" occupation probability r/ through 

When 77 < 77c, i?(7/) < 2"^?7'^+^ < r/. That is, iteration 
of the renormalization process will keep reducing the occu- 
pation probability further, as i?''+^(7/) < {rj) < . . . < rj. 
Better upper bounds on the renormalized occupation densities 
can be easily computed by solving the recurrence R^^^{rj) < 
2^[i?^(77)]''+^ with the initial condition R°{ri) = rj, which 
yields 

R''{t]) < 2^[i+(''+i)+---+('''+i)""'l?7(''+i)'' 

^ 2A[{k+ir-i]/k^{k+ir 

Note that the occurrence of errors on the coarse-grained in- 
teraction graph is also modeled by a site percolation process 
because the errors in any given procedure are assumed to oc- 
cur independently of the errors in all other procedures, and our 
renormalization transformation has been defined only in terms 
of the restriction of global error configurations to individual 
procedures. After r levels of concatenation, the effective error 
rate is exactly the site occupation probability on i?'"(Fc„), and 
is bounded from above by cj^^(ci77)*^^^^^'' , where ci, C2 are 
constants related to the particular quantum computation code 
used. The threshold condition on the error rate is therefore 
that 7/ < Cj^^. 

The necessary number of concatenation levels can now be 
determined in the usual way ||3l0]: if the unencoded circuit 
Co has N gates, then the encoded circuit Cr will have N r- 
procedures, so that if we want the net computation error to 
be less than some e > 0, the effective error rate must be at 
most e/N. This is guaranteed whenever R^{rf) < e/N, which 
yields r = polylog(A^/e) 1 10]. That is, the complexity of the 
fault-tolerant quantum circuit will be larger than that of the 
unencoded circuit by a polylogarithmic factor 

The above estimate of the threshold is quite crude. The 
exact value of 7/c can be determined by finding nontrivial fixed 
points of the renormalization transformation 77 Riv) (note 
that the trivial values 77 = 0, 1 are also fixed points of R). It 
is not hard to see that the graph of R has the "S-shape" shown 
in Fig. 12] The exact value rjc of the threshold error rate is 
precisely the nontrivial fixed point of R. One can also use the 
"staircase construction" pictured there to see how the effective 
error rate goes down with the number of concatenation levels. 

Further information on the renormalized error rate i?(77) can 
be obtained by means of general arguments of percolation the- 
ory. To this end let us observe the following properties of the 
event E: (i) its occurrence depends on the status (i.e., occu- 
pied or vacant) of a finite number a of vertices, and (ii) it is an 
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FIG. 2: Left: the "S-shaped" curve of the renormalization transfor- 
mation R. Right: the portion of the curve for < 77 < r^c and the 
"staircase construction" used to illustrate how the effective error rate 
is reduced under iteration of R. 



increasing event, i.e., stable under addition of more occupied 
vertices. Then the following differential inequality holds for 
< 77 < 1 liLlil: 



R{r^) (1 - R{t^)) 

T]{1 - 77) 



< R'iv) < 



77(1-77)' 



(1) 



where the prime denotes differentiation with respect to 77. 
It follows easily from the first inequality in Q that A = 
R'jrir) > 1; in fact, a more careful analysis shows that A > 1 
II13I1 . From the second inequality in Q we get the following 
relation holds between A, a, and rjc- 



A < 



(2) 



Consider now a quantum computer operating in the sub- 
threshold regime, but very close to the threshold, i.e., rjc — 5< 
rj < rjc for some small 5 > Q. The minimum number 
of concatenation levels necessary to bring the effective er- 
ror rate down to some suitable small value e (where we can 
assume that 5 <C 77c — e) is controlled by A (see the stair- 
case construction in Fig. |3. Linearizing R around rjc as 



R{r]) = ric + X{ri — rjc) + o{\r] — rjc]) and iterating this process 
r times, we see that we need to pick r such that rjc^ e < X^S, 
i.e., it suffices to take 



log[(77, - e)/S] 
log A 



(3) 



Solving (|3} for A, substituting into and using the fact that 

(77c — e)^/'' > (77c — e)^ for r > 1 and < e < 77^ < 1, we 
get 



Veil - Vc){Vc - < aS'^\ 



(4) 



This inequality can be regarded as a threshold-overhead trade- 
off inequality: it shows that fault-tolerant quantum circuits 
with low concatenation overhead must have correspondingly 
small threshold error rates, and conversely that fault-tolerant 
quantum circuits with large threshold rates have the distinct 
disadvantage of requiring high concatenation overhead. 

Let us close with a few remarks. First of all, we have so 
far paid no attention to the percolation threshold 77* on Tc 
and how it compares to 77c. We expect that rjc ^ 77^ on the 
grounds that the (connected) clusters of malfunctioning gates 
should be highly dilute in order for error correction to suc- 
ceed. Also, the methods used in this paper are applicable to 
noisy classical computers as well. (This is, ultimately, not 
very surprising in light of the fact that the local stochastic er- 
ror model is, essentially, classical in spirit.) It would be in- 
teresting and useful to extend the methods used in this paper 
(and percolation-theoretic methods in general) to less ideal- 
ized models of noise in quantum computers. A promising step 
in this direction would be an application of percolation tech- 
niques to problems involving graph-theoretic models of mul- 
tiparticle entanglement fl4] and quantum computation | fs"). 
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